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The simplest models of inflation predict small non-gaussianities and a featureless power spectrum. 
However, there exist a large number of well-motivated theoretical scenarios in which large non- 
gaussianties could be generated. In general, in these scenarios the primordial power spectrum 
will deviate from its standard power law shape. We study, in a model-independent manner, the 
constraints from future large scale structure surveys on the local non-gaussianity parameter /nl 
when the standard power law assumption for the primordial power spectrum is relaxed. If the 
analyses are restricted to the large scale-dependent bias induced in the linear matter power spectrum 
by non-gaussianites, the errors on the /nl parameter could be increased by 60% when exploiting 
data from the future DESI survey, if dealing with only one possible dark matter tracer. In the 
same context, a nontrivial bias |<5/nl| ~ 2.5 could be induced if future data are fitted to the 
wrong primordial power spectrum. Combining all the possible DESI objects slightly ameliorates the 
problem, as the forecasted errors on /nl would be degraded by 40% when relaxing the assumptions 
concerning the primordial power spectrum shape. Also the shift on the non-gaussianity parameter 
is reduced in this case, |<5/nl| ~ 1.6. The addition of Cosmic Microwave Background priors ensure 
robust future /nl bounds, as the forecasted errors obtained including these measurements are almost 
independent on the primordial power spectrum features, and |5/nl| ~ 0.2, close to the standard 
single-field slow-roll paradigm prediction. 

PACS numbers: 


I. INTRODUCTION 

Inflationary theories have been extremely successful in 
explaining the horizon problem and the generation of the 
primordial perturbations seeding the structures of our 
current universe um- The firm confirmation of these 
theories as the responsible ones for the universe we ob¬ 
serve today would come from the detection of a signal 
of primordial gravitational waves. A key observable to 
disentangle between different inflationary theories is the 
primordial power spectrum, i.e. the power spectrum of 
the initial curvature perturbations Pn(k). This power 
spectrum is usually taken to be a featureless primordial 
power spectrum (PPS), described by a simple power-law 
P-ji(k) oc k ris ~ 1 , with n s the scalar spectral index. How¬ 
ever, there exists a vast number of models in the liter¬ 
ature which may give rise to a non-standard PPS (see 
the recent review [T2]). That is the case of slow-roll in¬ 
duced by phase transitions in the early universe mama, 
by some inflation potentials HEHSZI, by resonant particle 
production [33IIT2L variation in the sound speed of adi¬ 
abatic modes [SI [44] or by trans-Planckian physics SU¬ 
ITS]- All these non-standard scenarios, as well as other 
non-canonical schemes [5tH457) , could lead to a PPS which 
may notably differ from the simple power-law parameter¬ 
ization. 

Another key observable to distinguish among the pos¬ 
sible inflationary models is the deviation from the pure 
Gaussian initial conditions. Non-gaussianities are usu¬ 
ally described by a single parameter, /nl- In the matter 


dominated universe, the gauge-invariant Bardeen poten¬ 
tial on large scales can be parametrized as [5M(jT| 

$ng = $ + /nl (d> 2 — (4> 2 )) , (1) 

where $ is a gaussian random field. The non-gaussianity 
parameter /nl is often considered to be a constant, yield¬ 
ing non-gaussianities of the local type. 

Traditionally, the standard observable to constrain 
non gaussianities is the Cosmic Microwave Background 
(CMB), through the three point correlation function, or 
bispectrum. As the odd power correlation functions van¬ 
ish for the case of Gaussian random variables, the bispec¬ 
trum provides the lowest order statistic to test any de¬ 
parture from gaussianity. The bispectrum is much richer 
than the power spectrum, as it depends on both the scale 
and the shape of the primordial perturbation spectra. 
The current bound from the complete Planck mission 
for the local non-gaussianity parameter is /nl = 0.8 ± 5 
(68% CL) (S2J. 

The large scale structures of the universe provide an 
independent tool to test primordial non-gaussianites, as 
shown in the pioneer works of Refs. [63] and [Mj . Dark 
matter halos will be affected by the presence of non- 
gaussianties, and a scale-dependent bias will characterise 
the non-gaussian signal at large scales IfioIfTl] . The 
tightest bounds on primordial non-Gaussianity using ex¬ 
clusively large scale structure data are those obtained 
from DR8 photometric data, see Ref. E2, which exploits 
800000 quasars and finds —49 < /nl < 31 (see also 
Ref- [73]). While current large scale structure constraints 
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are highly penalised due to their systematic uncertain¬ 
ties, it has been shown by a number of authors that the 
prospects from upcoming future large scale structure sur¬ 
veys can reach ct(/nl) < 1 i TTUTTHM] . 

Even if Eq. ([I]) is commonly used in the literature as¬ 
suming a scale independent parameter /nl , let us men¬ 
tion that some theoretical scenarios can give rise to a 
scale-dependent /nl [85H8B] . This scale dependence has 
already been studied in several works, see e.g. [59H9T] 
using large scale structure information, cluster number 
counts and/or CMB spectral information. The fore¬ 
cast on the errors on the non-gaussianity parameters 
are however known to be parametrization/model depen¬ 
dent [90]) 91]. The recent work of Ref. 91] focuses on 
the complementarity of the different cosmological probes, 
which could help enormously to determine the functional 
dependence of a scale-dependent non-gaussianity param¬ 
eter without having to assume a particular choice of such 
a scale-dependence. In particular, they make use of spec¬ 
tral distortions of the CMB background. In this work, 
we shall focus on the forecasts associated to future large 
scale structure probes only and we will restrict ourselves 
to a scale independent parameter /nl- However, when al¬ 
lowing for a non-standard primordial power spectrum as 
well, additional measurements of the CMB distortion pa¬ 
rameters could help in removing some of the degeneracies 
that appear between non-gaussianities and the parame¬ 
ters governing the primordial power spectrum parameter¬ 
ization. Furthermore, these degeneracies could copiously 
appear in the case of scale-dependent non-gaussinities. 

Despite the fact that the simplest models of inflation 
(i.e. single field, slow-rolling with a canonical kinetic 
term) predict small non-gaussianities, there are some the¬ 
oretical scenarios in which large non-gaussianties could 
be generated, see e.g Ref. (52] and references therein. 
The same deviations from the standard slow-roll infla¬ 
tion that give rise to non-gaussianities could also be 
a potential source for other features in the PPS m , 
which are absent in the simplest inflation models. Par¬ 
ticle production during inflation gives rise to both a 
non-canonical PPS and large non-gaussianities simulta¬ 
neously [32]. These two phenomena could also appear to¬ 
gether in single field models with non-standard inflation¬ 
ary potentials (2QU2I1 [24U32l[34l[36] , as well as in Brane 
Inflation [29] and multi-field inflationary models [55] , 
Other possibilities that will give rise to both a non¬ 
standard matter power spectrum and non-gaussianities 
include preheating scenarios 93 [94]. 

As nature could have chosen other inflationary scenario 
rather than the single field slow-roll paradigm, it is inter¬ 
esting to explore, in a model-independent way, how the 
forecasts for large scale structure surveys concerning fu¬ 
ture measurements of /nl are affected when the assump¬ 
tion of a standard PPS is relaxed. This has never been 
done before while forecasting errors on the /nl parameter 
and it is a mandatory calculation, because models which 
will produce non-gaussianities will likely give rise to a non 
standard PPS as well. Even if non-gaussianties and dis¬ 


tortions from the standard power-law PPS are expected 
to be governed by the same fundamental physics, (and 
therefore, related to each other), the underlying infla¬ 
tionary mechanism is unknown a priori. A conservative 
and general approach is therefore to treat these two phys¬ 
ical effects as independent and to be determined simul¬ 
taneously. This is the strategy we follow in this paper. 
The structure of this manuscript is as follows. We start 
describing the parameterization of the PSS used here in 
Sec. HU Section |III A| describes the scale-dependent halo 
bias in the matter power spectrum, while in Section [TTTB] 
we describe the methodology followed for our calcula¬ 
tions as well as the specifications of the future large scale 
structure survey illustrated here. We present our results 
in Sec. lIIICl and conclude in Sec. m 


II. PRIMORDIAL POWER SPECTRUM 

The simplest models of inflation predict a power-law 
form for the PPS of scalar and tensor perturbations. As 
previously stated, in principle, a different shape for the 
PPS (see Ref. [T2] and references therein), can be gen¬ 
erated by more complicated inflationary models (see e.g. 
Ref. [95] for some compilation). In order to explore the 
robustness of future forecasted errors from large scale 
structure surveys on the local non-gaussianity parame¬ 
ter /nLj we assume a non-parametric form for the PPS, 
following the prescription of Ref. [96], which is an ex¬ 
ample of a number of possible methods explored in the 
literature [9711125] . We describe the PPS of the scalar 
perturbations by means of a function to interpolate the 
PPS values in a series of nodes at fixed position. The 
function we exploit to interpolate is commonly named as 
a piecewise cubic Hermite interpolating polynomial , the 
PCHIP algorithm [126] . see the Appendix A of Ref. [96] 
for details concerning the version of the (PCHIP) algo¬ 
rithm m used in the following. Within this model, 
one only needs to provide the values of the PPS in a dis¬ 
crete number of nodes and to interpolate among them. 
As in previous work ,96], we define the PPS at twelve 
nodes, whose values of k are: 

k x = 5 • 10 -6 Mpc -1 , 
fc 2 = 10 -3 Mpc -1 , 

kj = fc 2 (fcii/fc 2 ) (j-2)/9 for j <£ [3,10], 
fcu = 0.35 Mpc -1 , 

k 12 = 10 Mpc -1 . (2) 

In the range (fc 2 , fen), that has been shown to be well con¬ 
strained by current cosmological data [124] . we choose 
equally spaced nodes (in logarithmic scale). The pur¬ 
pose of the first and the last nodes is to allow for a 
non-constant behaviour of the PPS outside the well- 
constrained range. The PCHIP PPS is given by 

P(fc) = P 0 x PCHIP(fc; P s ,i,..., P s ,i 2 ) , (3) 
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with P s j the value of the PPS at the node kj divided by 
Po = 2.2 • 10” 9 , according to the latest results from the 
Planck collaboration, see Ref. (128] . 


III. FORECASTS 
A. Non-gaussian halo bias 

Non-gaussianities as introduced in Eq. 0 induce a 
scale-dependent bias that affects the matter power spec¬ 
trum at large scales. This scale-dependent bias reads 
as [63] 165] 

8g — b ddm where b = be + Ab , (4) 

where <5 fl (<5dm) are galaxy (dark matter) overdensities, bo 
is the gaussian bias and A b reads as 

A 6 = 3/ N m-Wp||^, (5) 

with T(k) the linear transfer function. The growth factor 
D(a ) is defined as <Wi(a)/<Wi(a = 1) and 8 C refers to the 
linear overdensity for spherical collapse [129j . The power 
spectrum with non-gaussianties included is obtained us¬ 
ing 

P ns = P{b G + Ab + fnl) 2 , (6) 

where /i*, is the cosine of the angle between the line 
of sight and the wave vector k and / is defined as 
dlnddm/dlna. P is the dark matter power spectrum, 
whose k dependence is driven either by Eq. © or by the 
standard power-law matter power spectrum (with a given 
amplitude A s and slope n s ). 

In Fig. [IJ top panel, we illustrate the galaxy power 
spectrum in the absence of non-gaussianites (i.e. /nl = 
0) as well as for the /nl = 20 case. We also show, 
see the thin red dashed line, that using a PCHIP PPS 
with /nl = 0 it is possible to match the galaxy power 
spectrum obtained with a standard power-law PPS and 
/nl 7^ 0. The P s j values needed to obtain such an effect 
were taken to be within their currently 95% CL allowed 
regions [55]. Therefore, large degeneracies between the 
P s j nodes and /nl parameter are expected. Such a large 
value of /nl = 20, albeit allowed by the current large 
scale structure limits on local non-gaussianities, is much 
larger than the expected errors from the upcoming galaxy 
redshift surveys (see e.g. Refs. [751 T5U]). Therefore, we 
also illustrate in Fig. [I] bottom panel, the equivalent plot 
for /nl = 5. For this case, the values for the PPS nodes 
P s j required to match the predictions from a standard 
power law PPS lie within their 68% CL current allowed 
regions [96]. However, notice that the degeneracies are 
still present. We therefore expect that the forecasted er¬ 
rors on the /nl parameter are largely affected by the 
uncertainties on the precise PPS shape. 




k [/i/Mpc] 


FIG. I: The upper panel depicts the galaxy power spectrum 
for the standard PPS power law case, for /nl = 0 (black solid 
curve) and /nl = 20 (blue dotted curve), together with a 
PCHIP PPS case (red dashed lines) for /nl = 0. The values 
of the PCHIP PPS nodes are chosen accordingly to match the 
predictions of the /nl = 20 case. The lower panel shows the 
equivalent but for /nl = 5. We have also changed accordingly 
the value of the PCHIP PPS nodes. We also show with ki for 
i = 4,.., 11 the k position of five of the nodes considered in our 
analysis (i = 5,.., 9), plus A? 4 ,io,ii that lie outside the k range 
probed by the DESI experiment. The galaxy power spectra 
are obtained for z = 0.57, \/j,k\ = 1 and assuming a constant 
gaussian bias &g. 


B. Methodology 

We focus here on the future spectroscopic galaxy sur¬ 
vey DESI (Dark Energy Instrument) experiment [ 1311 . 
Although multi-band, full-sky imaging surveys have 
been shown to be the optimal setups to constrain non- 
gaussianities via large scale structure measurements cm 
E5, the purpose of the current paper is to explore the de¬ 
generacies with the PPS parameterization rather than to 
optimise the /nl sensitivity. Therefore, we restrict our¬ 
selves here to the DESI galaxy redshift survey (similar 
results could be obtained with the ESA Euclid instru¬ 
ment) . 

In order to compute the expected errors on the lo¬ 
cal non-gaussianity parameter, we follow here the usual 
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Fisher matrix approach, whose elements, as long as the 
posterior distribution for the parameters can be well ap¬ 
proximated by a Gaussian function, read as |132H134| 

F a ^ = ^Tr[C- 1 C <a C- 1 C^\ , (7) 


with C = S+N the total covariance. The covariance ma¬ 
trix contains both the signal S and the noise N terms, 
and C\ a refer to its derivatives with respect to the cos¬ 
mological parameter p a in the context of the underlying 
fiducial cosmology. The 68 % CL marginalized error on a 
given parameter p a is u{p a ) = \/ (F _ 1 )a a , with F _1 the 
inverse of the Fisher matrix. In the following, in order to 
highlight the differences in the error on the /nl param¬ 
eter arising from different PPS choices, we only consider 
information concerning non-gaussianites from large scale 
structure data, neglecting the information that could be 
added from CMB bispectrum measurements. 

Our large scale structure Fisher matrix reads as [135) 


771LSS _ 

*at(3 ~ 


L 


d In P ng {k) d In P ng {k) 


Kff(fc); 


dk 


1 rk„ 


1-1 Jk„ 


dp a dpp 2(2tt) 3 

1 d\nP ng (k,p k ) 3In P ng (k, p k ) 


( 8 ) 


dp 0 


dpp 


V e s(k, Pk) 


2nk 2 dkdpk 

2(2t r) 3 


where T 4 ff is the effective volume of the survey: 


Veff{k : Pk) 


nP n g(k,fj k ) ] 2 
nP ng (k, p,k) + 1. 


(9) 


where P ng is the power spectrum with non-gaussianities 
included (see Eq. ([ 6 ])) and n refers to the galaxy number 
density per redshift bin. We assume fc max = 0.1/i/Mpc 
and fc m i n is chosen to be equal to 2tt/V 1 ^ 3 , where V rep¬ 
resents the volume of the redshift bin. The DESI survey 
is expected to cover 14000 deg 2 of the sky in the redshift 
range 0.15 < 2 < 1.85, divided in redshift bins of width 
A z = 0.1. We follow Ref. m for the number densities 
n(z) and biases bc(z) associated to the three types of 
DESI tracers: Luminous Red Galaxies (LRGs), Emission 
Line Galaxies (ELGs) and high redshift quasars (QSOs). 
We include the redshift dependence of the (fiducial) bias 
bo in Eq. Q as follows: bc{z)D{z) = 0.84,1.7,1.2 for 
ELG, LRG and QSO’s respectively, where D{z) is the 
growth factor as in Eq. (|5j). In order to combine the three 
different Fishers matrices from the three DESI tracers 
(LRGs, ELGs and QSOs), we follow the multi-tracer for¬ 
malism developed in Ref. D323, where the authors present 
a generic expression for the Fisher information matrix of 
surveys with any number of tracers. The multi-tracer 
technique provides constraints that can surpass those set 
by cosmic variance, due to the differences in the cluster¬ 
ing of the possible tracers of large scale structure. 

Also, let us remind that the observed size of an object 
or a feature at a given redshift z are obtained in terms of 


redshift and angular quantities A z and A 9. These two 
quantities are related to the comoving dimensions rn and 
r _l along and across the line of sight through the angular 
diameter distance Da{z) and the Hubble rate H(z). The 
same applies to the Fourier transform associated vari¬ 
ables (we will refer to these as fen and k± for the dual 
coordinates of ry and rj_). Therefore, when one aims to 
reconstruct the measurements of galaxy redshifts and po¬ 
sitions in some reference cosmological model which differ 
from a given fiducial cosmology, one has to account for 
geometrical effects in the following way [135 ] : 


Po bs (k7 f ,kl ef ) 


D a(z) \ 2 ef 

D a (z) 2 


H(z) 

H{z)\ re f 


Pf id{k\\ , k_ L) 


. (10) 

where the ref sub/superscript denote quantities in the 
reference cosmological model*. We properly account for 
such effects in our Fisher matrix forecasts when taking 
numerical derivatives of the galaxy power spectrum with 
respect to the cosmological parameters at given values of 
|k| and p, k (or equivalently, of fcy and k±). 

In addition to Fisher matrix forecasts, we will also com¬ 
pute the expected shift in the /nl parameter if the P s j 
PCHIP parameters (with j = 5, ..,9) are (incorrectly) set 
to values different from their fiducial ones. For that pur¬ 
pose, we use the method developed by the authors of 
Ref. [1138 :. The main idea is as follows: if the future DESI 
data are fitted assuming a cosmological model with fixed 
values of P s ,j' and therefore characterised by n' = 5 pa¬ 
rameters M! = {flbh 2 , fl c h 2 , h, /nl, w}, but the true 
underlying cosmology is a model with different values 
of P s j and therefore characterized by n = 10 parameters 
M = {fl b h 2 , Q c h 2 , h, /nl, w, P s j} (with j = 5, ..,9), the 
inferred values of the n! = 5 parameters will be shifted 
from their true values to compensate for the fact that 
the model used to fit the data is wrong. Assuming that 
the likelihood is gaussian, the shifts in the n' parameters 
read as m 


8e' a = -{F' 1 ) a pG K 5ip< : a,0=l...nf, 

C = n' + 1... n , (11) 


where F' is the Fisher matrix for the n' parameters model 
(with the P s j fixed) and G denotes the Fisher matrix 
for the n parameters model (including the n! previous 
parameters plus the PCHIP P S J parameters). 

In the following, unless otherwise stated, we shall 
adopt the best-fit values from the complete Planck mis¬ 
sion [1281 . which, in the standard power law PPS, corre¬ 
sponds to A s = 2.2-10 -9 and n s = 0.965 at k plvot = 0.05. 
Within the PCHIP parameterization, the best-fit values 
used for the nodes considered in the numerical analy¬ 
ses below are: P S} 5 = 1.07099, P S fi = 1.04687, P s> 7 = 


* fc ll = k \\ ef D A(z)\ re f/D A (z) and k ± = k r f f H(z)/H(z)\ TB f . 
t Fixing the values of P s j corresponds to fix both n s and A s to 
their best-fit values according to the normalization used here. 
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1.02329, P s , 8 = 1.00024 and P s , 9 = 0.97771. These 
values are obtained calculating the value of the best-fit 
power-law power spectrum, given by Planck 2015 best-fit 
values for A s and n s as mentioned above, at the positions 
of the nodes kg to kg. The remaining nodes are outside 
the k range expected to be covered by the DESI survey, 
given the values of fc max and fc m ; n considered here. 


C. Results 

In the following, we shall present the results arising 
from our Fisher matrix calculations, for the two fiducial 
cosmologies explored here: one in which the PPS is de¬ 
scribed by its standard power-law form, and a second 
one in which the PPS is described by the PCHIP param¬ 
eterization. The parameters describing the model with 
a power-law PPS are the baryon and cold dark matter 
physical energy densities, f Igh 2 and fi c ft 2 , the Hubble 
parameter h (with Hg, the Hubble constant, defined as 
100ft. km/s/Mpc), the scalar spectral index n s , the am¬ 
plitude of the PPS A s , and the equation of state of the 
dark energy component w. The PCHIP PPS case is also 
described by H&ft 2 , f2 c ft 2 , ft, w plus five nodes P s j (with 
j ranging from 5 to 9) describing the PCHIP PPS. Non- 
gaussianites of the local type are implemented in both 
fiducial cosmologies via the /nl parameter. All the re¬ 
sults described below (unless otherwise stated) refer to 
the analysis of the three DESI tracers (ELGs, LRGs and 
QSOs), i.e. they have been obtained exploiting exclu¬ 
sively the scale-dependent biases imprinted in the power 
spectra of these three types of tracers. 

Table H<0> shows the lcr marginalized errors for the 
case of a standard (PCHIP) PPS, for a fiducial value 
/nl = 20 for each of the DESI tracers as well as the er¬ 
ror from the combination of all of them, using the multi¬ 
tracer technique. Even if such a value of the /nl param¬ 
eter (/nl = 20 ) is larger than the expected sensitivity 
from future probes, it is still allowed by current large 
scale structure bounds on primordial non-gaussianities. 
Notice that, for the standard power law PPS, the ex¬ 
pected error on /nl is 19.9, 10.1 and 8.56 for LRGs, 
ELGs and QSOs respectively, while for the case of the 
PCHIP parameterization, one obtains ct(/nl) = 32.2, 13.3 
and 12.6 respectively. Therefore, there is a large increase 
in the error on the non-gaussianity parameter, which can 
reach the 60% level. Concerning the remaining cosmolog¬ 
ical parameters, they are barely affected. In some cases, 
their error is even smaller than in the standard power- 
law scenario. This is indeed the case of the equation 
of state parameter w, or Hf,ft 2 and f l c h 2 (the errors on 
the latter two parameters are smaller than in the stan¬ 
dard PPS approach only when exploiting either ELGs or 
QSOs tracers). The combination of the data from the 
three tracers exploiting the multi-tracer technique allevi¬ 
ates the problem with the error on /nl> as the increase 
in the value of ct(/nl) when relaxing the assumption of a 
simple power-law PPS is around 40%, rather than 60%. 


The reason for this generic increase in the error of /nl 
is due to the the large degeneracies between the non- 
gaussianity /nl parameter and the P S j nodes, which get 
reduced when combining the tracers. The top and bot¬ 
tom panels of Fig. [TT] illustrate the large degeneracies be¬ 
tween the non-gaussianity /nl parameter, for the fiducial 
value /nl = 20 and two of the P s j nodes, P Sj5 and P Sj 9 . 
We only show here these two nodes, but similar degen¬ 
eracies are obtained for the remaining nodes. 


This degeneracy problem could a priori be solved in 
two ways, either exploiting smaller scales in the observed 
galaxy or quasar power spectra, or using CMB priors. 
In practice, going to the mildly non-linear regime would 
require new additional P s j nodes and new degenera¬ 
cies between these additional P s j nodes and the non- 
gaussianity parameter /nl will appear. We have numer¬ 
ically checked that such a possibility does indeed not 
solve the problem. Furthermore, a non-linear descrip¬ 
tion of the matter power spectrum will depend on addi¬ 
tional parameters, enlarging the number of degeneracies. 
In contrast, the CMB priors on both the PPS param¬ 
eters as well as on the dark matter and baryon mass- 
energy densities help enormously in solving the problem 
of the large degeneracies between the PPS parameteriza¬ 
tion and non-gaussianities. Tables m and |IV| show the 
equivalent of [I] and [TT] but including CMB priors from 
the Planck mission 2013 data [139] . Notice that the im¬ 
pact of the Planck priors is largely more significant in the 
PCHIP parameterization case: the /nl errors arising from 
the three different dark matter tracers when the CMB 
information is included are smaller in the PCHIP PSS de¬ 
scription than in the standard power-law PSS modeling. 
When the multi-tracer technique is applied, the overall 
errors after considering Planck 2013 CMB constraints are 
very similar regardless on the PPS description and close 
to <t(/nl) - 5. 


Table [V (VII shows the lcr marginalized errors for the 
case of a standard (PCHIP) PPS, for another possible non- 
gaussianity parameter fiducial value, /nl = 5, from each 
of the DESI tracers, as well as the error arising from 
the combination of all of them using the multi-tracer 
technique. As in the case of /nl = 20, the error on 
the non-gaussianity parameter is increased, reaching in 
some cases a 60% increment. The results are very similar 
to those obtained and illustrated before for larger non- 
gaussianites. The errors on the other cosmological pa¬ 
rameters remain unaffected under the choice of the PPS 
parameterization. The dark energy equation of state pa¬ 
rameter is extracted with a smaller error in the PCHIP 
PPS case, and also f Igh 2 and H c ft 2 are determined with a 
smaller error in that case while dealing with either ELGs 
or QSOs tracers. The multi-tracer technique provides a 
reduction on the /nl error similar to that obtained in 
the /nl = 20 case. The top and bottom right panels of 
Fig- E| illustrate the large degeneracies between the non- 
gaussianity /nl parameter and two of the P s j nodes, 
P Sj 5 and P Sj9 for the fiducial value /nl = 5. Notice that 
the degeneracy pattern appears to be independent of the 
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fiducial 

LRG 

ELG 

QSO 

All 


0.02267 

4.78 • 10~ 3 

4.86 • 10 -3 

5.11 ■ 10~ 3 

2.38- 10“ 3 

Q c h 2 

0.1131 

1.75- 10" 2 

1.65 ■ 10” 2 

1.51 ■ 10~ 2 

7.70- 10“ 3 

h 

0.705 

5.02 • 10“ 2 

5.01 • 10' 2 

4.69 • 10' 2 

2.42 • 10~ 2 

n s 

0.96 

5.68 • 10~ 2 

4.28 • 10~ 2 

4.12 • 10~ 2 

1.96- 10" 2 

A s 

2.2- 10" 9 

0.341 

0.331 

0.302 

0.156 

/nl 

20 

19.9 

10.1 

8.56 

4.79 

w 

-1 

5.38- 10~ 2 

4.09 ■ 10“ 2 

6.18- 10~ 2 

2.36- 10 -2 


TABLE I: Marginalized l-cr constraints on the parameters associated to the standard PPS assuming a fiducial value /nl = 20. 
The error on the amplitude of the power spectrum is evaluated on A a /( 2.2 • 10 -9 ). 



fiducial 

LRG 

ELG 

QSO 

All 


0.02267 

7.85- 10~ 3 

3.65 • 10~ 3 

4.70 • 10 -3 

2.30 ■ 10~ 3 

Q c h 2 

0.1131 

2.30- 10~ 2 

1.11 ■ 10~ 2 

1.41 ■ 10 -2 

6.36- 10~ 3 

h 

0.705 

7.67- 10~ 2 

3.59- 10" 2 

4.62 • 10 -2 

2.12 • 10~ 2 

Ps,5 

1.07099 

0.340 

0.169 

0.212 

0.111 

Ps, 6 

1.04687 

0.419 

0.198 

0.254 

0.119 

Ps,7 

1.02329 

0.451 

0.216 

0.276 

0.125 

Ps, 8 

1.00024 

0.479 

0.229 

0.293 

0.132 

Ps, 9 

0.97771 

0.482 

0.234 

0.298 

0.134 

/nl 

20 

32.2 

13.3 

12.6 

6.43 

w 

-1 

4.03- 10~ 2 

2.80 • 10~ 2 

4.45 • 10 -2 

2.45 ■ 10~ 2 


TABLE II: Marginalized l-cr constraints on the parameters associated to the non-standard PPS assuming /nl = 20. 


value of /nl- The addition of the CMB priors brings the 
errors on all the cosmological parameters (/nl included) 
to the same values in both PPS parameterizations (stan¬ 
dard power-law and PCHIP PPS prescriptions), as shown 
in Tabs. lyTH and IVTTn 


We now perform an additional forecast. We focus here 
on the shift induced in the local non-gaussianity param¬ 
eter /nl, which we set to zero in the two cosmologies M 
and AT. For the purpose of this analysis, we fix all of the 
P s j to their best-fit values according to the Planck 2013 
results, for the case of the AT cosmology. A shift in /nl 
is expected to compensate for the fact that the P S} j PCHIP 
nodes are additional parameters in M, while not being 
considered as free parameters in the AT analysis. There¬ 
fore, we displace the P s j parameters (with j = 5,..,9) 
from their fixed fiducial values in AT, i.e. we are adding 
them as additional parameters in the cosmological model 
(i.e. to be determined by data). Referring to the no¬ 
tations of Eq. |TT|, using a shift =0.1, which is 

smaller than theTr expected errors (see Tabs. Pand[Vl|, 
we obtain that the corresponding shift in the /nl parame¬ 
ter is Sdf KL ~ 2.5, regardless of the exploited dark matter 
tracer. This is a quite large displacement of the local non- 
gaussianity parameter which will induce a non-negligible 
bias in reconstructing the inflationary mechanism. While 
the remaining cosmological parameters are also slightly 
displaced with respect to their fiducial values, their shifts 
will not induce a misinterpretation of the underlying true 
cosmology. The non-gaussianity shift could be a 

potential problem when extracting the (true) value of 
the /nl parameter not only for the DESI survey, but 
also for future experiments with improved sensitivities 
to non-gaussianities, such as SPHEREx m ■ The com¬ 


bination of all the three possible DESI tracers leads to 
a smaller shift in the /nl parameter {Sdf NIj ~ 1.6). If 
CMB priors are applied, the shift is considerably reduced, 
<50 / nl ~ 0 .2, which is close to the expectations for non- 
gaussianities in the most economical inflationary models, 
i.e. within single field slow-roll inflation [9211140] . 


IV. CONCLUSIONS 

While the simplest inflationary picture describes the 
power spectrum of the initial curvature perturbations 
Pp(k) by a simple power-law without features, there ex¬ 
ists a large number of well-motivated inflation models 
that could give rise to a non-standard PPS. The major¬ 
ity of these models will also generate non-gaussianities. 
The large scale structure of the universe provides, to¬ 
gether with the CMB bispectrum, a tool to test primor¬ 
dial non-gaussianites. Plenty of work has been devoted in 
the literature to forecast the expectations from upcom¬ 
ing galaxy surveys, such as the Dark Energy Instrument 
(DESI) experiment. The forecasted errors and bounds on 
the non-gaussianity local parameter /nl are however usu¬ 
ally derived under the assumption of a standard power- 
law PPS. Here we relax such an assumption and compute 
the expected sensitivity to /nl from the DESI experi¬ 
ment assuming that both the precise shape of the primor¬ 
dial power spectrum and the non-gaussianity parameter 
need to be extracted simultaneously. If the analysis is 
restricted to large scale structure data, the standard er¬ 
rors computed assuming a featureless power spectrum are 
enlarged by 60% within the PCHIP PPS parameterization 
explored here and when treating each of the possible dark 
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fiducial 

LRG ELG QSO All 


0.02267 

2.67 ■ 10~ 4 2.63 ■ 10 -4 2.66 • 10” 4 2.59 • 10“ 4 

Q c h 1 2 

0.1131 

1.64 • 10 -3 1.44 • 10“ 3 1.52 • 10~ 3 1.24 • 10 -3 

h 

0.705 

6.66 • 10" 3 5.24 • 10 -3 5.86 • 10~ 3 4.12 • 10“ 3 

n s 

0.96 

6.72 • 10 -2 6.41 • 10~ 2 6.53 • 10~ 2 5.84 • 10 -3 

A s 

2.2- 10" 9 

3.87 ■ 10“ 2 3.28 • 10 -2 3.51 ■ 10“ 2 2.71 • 10 -2 

/nl 

20 

17.4 9.14 7.58 4.56 

w 

-1 

4.51 ■ 10 -2 3.36 • 10“ 2 5.44 • 10~ 2 2.17 • 10 -2 


TABLE III: As Tab. [I] but including CMB priors, see the text for details. 



fiducial 

LRG 

ELG 

QSO 

all 

Q.bh 2 

0.02267 

3.92- 

10~ 4 

3.79- 

10“ 4 

3.87 

• 10" 4 

3.74 

10“ 4 

Q c h 2 

0.1131 

1.36- 

10” 3 

1.10- 

10“ 3 

1.18 

• 10 -3 

1.04 

10 -3 

h 

0.705 

4.13- 

10~ 3 

3.14- 

10“ 3 

3.62 

• 10 -3 

2.93 

10 -3 

Ps, 5 

1.07099 

2.98- 

10” 2 

2.69- 

10 -2 

2.77 

• 10 -2 

2.60 

10 -2 

Ps, 6 

1.04687 

2.89- 

10~ 2 

2.10 • 

10“ 2 

2.32 

• 10 -2 

1.99 

10 -2 

Ps, 7 

1.02329 

2.00- 

10~ 2 

1.73- 

10 -2 

1.84 

• 10 -2 

1.69 

10 -2 

Ps, 8 

1.00024 

1.92 • 

10~ 2 

1.76- 

10“ 2 

1.86 

• 10' 2 

1.73 

10' 2 

Ps, 9 

0.97771 

2.59- 

10~ 2 

2.31 ■ 

10 -2 

2.42 

• 10 -2 

2.22 

10 -2 

/nl 

20 

13.0 

6.85 

5.64 

4 

75 

w 

-1 

3.24- 

10~ 2 

2.46 • 

10“ 2 

4.0- 

10~ 2 

2.28 

10 -2 


TABLE IV: As Table |H] but including CMB priors. 


matter tracers individually. Another potential problem 
in future galaxy surveys could be induced by the (wrong) 
assumption of a featureless PSS (while nature could have 
chosen a more complicated inflationary mechanism lead¬ 
ing to a non-trivial PPS). If future data is fitted to the 
wrong PPS cosmology, a shift in |<50/ NL | ~ 2.5 would be 
inferred (for k max = O.lh/Mpc) even if the true cosmol¬ 
ogy has /nl = 0. The multi-tracer technique helps in 
alleviating the former two problems. After combining all 
the DESI possible tracers, the forecasted errors on /nl 
will be degraded by 40% (when compared to the value 
obtained within the standard power-law PPS model) and 
the resulting shift will be reduced to |<$0/ NL | ~ 1.6. The 
addition of Cosmic Microwave Background priors from 
the Planck 2013 data on the PPS parameters and on the 
dark matter and baryon mass-energy densities lead to 
a /nl error which is independent of the PPS parame¬ 
terization used in the analysis. After considering CMB 
priors, the value of the shift |<50/ NL | is 0.2, which is of 
the order of standard predictions for single-held slow-roll 
inflation iTl. 1140] . 
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fiducial 

LRG 
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QSO 

All 


0.02267 

4.78 • 10~ 3 
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5.18 • 10~ 3 
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Q c h 2 

0.1131 
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fiducial 

LRG 

ELG 

QSO 

All 


0.02267 

2.67- 

10~ 4 

2.63 ■ 10 -4 

2.67- 10~ 4 

2.59- 10“ 4 

f l c h 2 
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1.64 • 

10 -3 
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TABLE VII: As Tab. [V] but including CMB priors, see the text for details. 



fiducial 

LRG 

ELG 

QSO 

all 

flbh 2 

0.02267 

3.92- 

10~ 4 

3.79- 

10" 4 

3.86 

10" 4 

3.75- 10~ 4 

Q c h 2 

0.1131 

1.36- 

10~ 3 

1.10- 

10“ 3 

1.18 

10" 3 

1.04- 10~ 3 

h 

0.705 

4.10- 

10~ 3 

3.13- 

10“ 3 

3.59 

10" 3 

2.92 • 10~ 3 

Ps, 5 

1.07099 

2.98- 

10~ 2 

2.68- 

10“ 2 

2.77 

10" 2 

2.60 ■ 10~ 2 

Ps,6 

1.04687 

2.89- 

10“ 2 

2.11 ■ 

10“ 2 

2.33 

10" 2 

2.0 ■ 10“ 2 

Ps, 7 

1.02329 

2.00- 

10~ 2 

1.73- 

10“ 2 

1.84 

10" 2 

1.69- 10~ 2 

Ps, 8 

1.00024 

1.92 • 

10“ 2 

1.76- 

10“ 2 

1.86 

10" 2 

1.73 • 10~ 2 

Ps. 9 

0.97771 

2.50- 

10“ 2 

2.31 ■ 

10“ 2 

2.43 

10" 2 

2.22 ■ 10~ 2 

/nl 

5 

12.4 

6.42 

5 

23 

4.46 

w 

-1 

3.23- 

10~ 2 

2.46 • 

10“ 2 

3.99 

10 -2 

2.27- 10~ 2 


TABLE VIII: As Table [VT| but including CMB priors. 
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FIG. II: The upper left (right) panel shows the /nl — Ps, 5 degeneracy, for a fiducial cosmology with /nl = 20 (/nl = 5), 
assuming fc max = 0.1/i/Mpc. We show the 1 — a marginalized contours associated to the LRGs (in dashed blue lines), ELGs 
(in dot-dashed green lines), QSOs (in dotted cyan lines) and multi-tracer (in solid red) Fisher matrix analyses. The bottom 
panels shows the analogous but in the (/nl — Ps , 9) plane. 
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